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Abstract 


The  objective  of  the  study  conducted  under  this  grant  is  to  develop  effective  computational 
strategies  for:  a)  assessing  the  sensitivity  of  the  dynamic  failure  and  damage  of  metallic 
structures  to  variations  in  both  the  microstructure  and  macrostructure  material  parameters;  and  b) 
incorporating  the  uncertainties  in  material  parameters  in  the  prediction  of  response  and  damage 
of  these  structures.  The  sensitivity  coefficients,  which  are  the  derivatives  of  the  various 
structural  response  quantities  with  respect  to  the  material  parameters,  are  evaluated  by  using  the 
direct  differentiation  method. 

A  two-phase  computational  procedure  is  applied  to  the  uncertainty  analysis  of  dynamic  failure  in 
metallic  structures.  In  the  first  phase  the  sensitivity  analysis  is  used  to  identify  the  major 
material  parameters,  i.e.,  the  parameters  which  have  the  most  impact  on  the  response  and  failure. 
In  the  second  phase,  the  major  parameters  are  treated  as  fuzzy  parameters,  and  the  fuzzy  set 
techniques  are  used  to  provide  information  about  the  range  and  variation  of  possible  response 
quantities  and  damage  associated  with  pre-selected  ranges  and  variations  of  the  major  parameter. 

The  computational  procedure  developed  under  this  grant  can  be  used  to  generate  a  toolkit  of  the 
uncertainty  analysis  procedure  that  allows  commercial  deterministic  analysis  codes  to  perform 
uncertainty  analysis.  This  is  accomplished  by  modifying  the  input  (to  include  the  pre-selected 
variations  in  the  major  parameters),  and  adding  a  postprocessor  to  generate  the  range  of  variation 
of  different  response  quantities. 

A  number  of  thermo-viscoplastic  material  models  will  be  developed  to  study  the  effect  of  rate- 
sensitivity  on  welding  residual  stresses  in  stainless  steel  AL6XN.  Also,  material  models  are 
developed  which  reduce  to  viscous  models  above  the  solidus  temperature  Ts  and  to  rate 
independent  plastic  models  below  0.4  Ts.  Between  those  two  temperature  ranges,  the  models 
describe  rate-dependent  plastic  behavior.  Attempt  was  made  to  estimate  the  values  of  the 
material  parameters,  over  wide  ranges  of  temperatures  and  strain  rates,  by  working  closely  with 
researchers  performing  experiments.  A  number  of  arc  and  laser  welding  techniques  for  residual 
stress  mitigation  are  considered.  The  effect  of  variation  in  material  parameters,  which  have  the 
most  impact  on  the  response  and  damage,  are  studied  by  using  fuzzy-set  techniques. 
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1.  Background  and  Brief  Overview  of  Previous  Work 


The  formation  of  residual  stresses  is  inevitable  in  any  welding  operation.  When  combined  with 
the  externally  applied  stress  fields,  the  residual  stresses  can  result  in  significant  differences  in  the 
performance  of  welded  structures  such  as  dimensional  stability,  fatigue  strength,  and  resistance 
to  environmentally  induced  cracking  [Ref.  1],  Furthermore,  it  has  been  shown  [e.g.,  Refs.  2,  3] 
that  the  residual  stresses  can  significantly  affect  the  fracture  behavior  of  welded  components. 

The  formation  of  welding  residual  stress  in  engineering  materials  is  a  complex  thermal- 
mechanical-metallurgical  process.  A  welding  operation  introduces  localized  heating  and  cooling 
around  the  joint  region.  Material  expands  on  heating  and  contracts  on  cooling.  The  localized 
material  expansion  and  contraction  during  welding  create  a  transient,  non-uniform  deformation 
process  which  results  in  the  formation  of  residual  stresses  due  to  the  uneven  distribution  of  the 
inelastic  strains.  For  many  engineering  alloys,  the  welding  thermal  cycle  would  induce 
metallurgical  changes  in  the  heat  affected  zone  (HAZ)  and  the  weld  metal.  These 
microstructural  changes  could  have  profound  influence  on  the  distribution  of  residual  stresses  in 
the  weldment,  due  to  the  formation  of  complex  material  property  gradients,  the  transformation- 
induced  plasticity,  and  the  specific  volume  changes  of  the  microstructure  constituents  associated 
with  the  phase  transformations  in  the  HAZ  and  weld  metal.  The  thermal-mechanical- 
metallurgical  interactions  of  a  material  during  welding  and  their  effects  on  residual  stresses  are 
illustrated  in  Figure  1. 

As  shown  in  Figure  1,  the  material  to  be  welded  and  the  welding  process  conditions  are  two 
primary  factors  governing  the  formation  of  residual  stresses  of  a  weldment. 

The  welding  induced  residual  stresses  strongly  depend  on  the  material  to  be  welded.  Both  the 
thermophysical  properties  and  the  mechanical  properties  of  a  material  would  influence  the 
residual  stresses  in  a  welded  structure.  The  thermophysical  properties  (thermal  conductivity, 
specific  heat,  latent  heat  of  fusion,  and  density)  influence  the  temperature  field  surrounding  the 
weld  pool,  thereby  influencing  the  thermal  expansion  and  contraction  behaviors  of  the  material 
during  welding.  The  mechanical  properties  (elastic  modulus,  yield  strength,  flow  stress,  and 
thermal  expansion  coefficient)  contribute  directly  to  the  stress  evolution  during  a  welding  cycle. 
The  microstructural  and  associate  metallurgical  changes  contribute  property  changes  in  the  weld 
region  are  highly  dependent  on  the  type  of  alloy  system  and  the  composition  of  alloys. 

The  weld  residual  stresses  also  strongly  depend  on  welding  process  conditions.  Fundamentally, 
it  is  the  thermal  energy  from  a  welding  heat  source  that  drives  the  heat  flow  in  the  material  being 
welded.  The  spatial  and  temporal  distributions  of  the  welding  temperature  field  is  directly 
related  to  the  welding  process  conditions  such  as  welding  speed,  and  the  intensity  and 
distribution  of  the  heat  input.  Changes  in  the  temperature  field  not  only  alter  the  thermal 
expansion  and  contraction  process,  but  also  the  phase  transformation  behaviors  of  the  material 
during  welding;  both  would  influence  the  residual  stresses  and  distortions  of  welded  structures. 
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Welding  Process 
&  Parameters 


Strength,  Soundness  of  Weldment) 

Figure  1  Primary  factors  and  the  thermal-mechanical-metallurgical  interactions  during 
welding  that  contribute  to  the  formation  of  weld  residual  stresses  and 
performance  of  welded  structures. 


Understanding  the  effects  of  welding  process  conditions  on  weld  residual  stresses  is  not  only  a 
fundamental  subject  of  welding  science,  but  also  has  broad  practical  implications.  For  many 
Navy  applications,  the  selection  of  materials  and  the  design  of  a  welded  structure  are  often 
determined  based  on  the  overall  performance  requirements  of  the  structure;  they  are  seldom 
changeable  because  of  welding  related  considerations.  Thus,  from  a  practical  point  of  view, 
effective  control  and  mitigation  of  weld  residual  stresses  and  distortions  would  mainly  rely  on 
proper  selection  of  welding  process  conditions.  Over  the  years,  the  welding  community  has 
considered  many  welding  and  post-welding  techniques  to  manage  the  residual  stresses  and 
distortions  [Ref.  1].  However,  much  of  the  development  has  been  based  on  the  time-consuming 
and  costly  trial-and-error  experimental  testing,  and  often  with  limited  success. 

In  recent  years,  there  have  been  attempts  to  “pro-actively”  manage  the  thermomechanical 
deformation  process  during  welding  to  control  the  welding  residual  stresses  and  distortions. 
Among  the  proposed  approaches  are  the  so-called  low-stress  and  no-distortion  welding  for  thin 
structure  panels  [Ref.  4],  the  thermal  tensioning  welding  (a  derivative  of  the  no-distortion  and 
low-stress  technique)  for  shipyard  panels  [Ref.  5],  the  active  quenching  welding  [Ref.  6],  and  the 
use  of  low-temperature  phase  transformation  filler  metal  in  welding  [Ref.  7].  Although  limited 
success  have  been  reported  of  these  techniques,  they  are  still  in  the  early  stages  of  research  and 
development  -  the  process  window  in  which  each  of  the  aforementioned  techniques  operates  is 
generally  not  well  defined  in  terms  of  the  key  process  conditions  and  structural  geometry. 
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Clearly,  understanding  the  sensitivity  of  thermal-mechanical-metallurgical  responses  of  material 
to  variations  in  welding  process  conditions  would  be  a  critical  aspect  toward  widespread 
applications  of  the  above  welding  techniques  to  mitigate  the  welding  residual  stresses  and 
minimize  their  detrimental  effects.  Built  upon  the  progress  of  the  current  ONR  grant,  the  focus 
of  future  research  should  be  on  investigating  the  effects  of  the  welding  process  conditions. 

The  study  under  this  ONR  grant,  focused  on  the  development  of  an  integrated  computational 
welding  simulation  environment  that  can  be  used  to  systematically  investigate  and  interrogate  the 
sensitivity  and  variability  of  the  thermal-mechanical-metallurgical  responses  of  welded 
components  to  various  material,  geometric  and  welding  conditions.  The  computational 
simulation  environment  consists  of  (1)  advanced  welding  simulation  procedures,  (2)  effective 
sensitivity  and  fuzzy-set  analysis  methodology,  and  (3)  advanced  computer  visualization 
technology  to  investigate  and  interrogate  the  sensitivity  and  variability  of  the  response  to 
variations  in  material  parameters. 

The  advanced  welding  simulation  procedure  integrated  the  metallurgical  and  microstructural 
computation  algorithms  into  finite  element  based  thermo-mechanical  code.  This  allows  for 
fundamental  investigation  of  the  effects  of  material  properties  variations/gradients  in  the  weld 
joint  region  on  the  performance  of  welded  metallic  structures. 

The  framework  of  the  computational  welding  simulation  environment  was  developed  and 
applied  to  investigate  the  effect  of  uncertainties  in  material  properties. 


2.  Mathematical  Formulation 


This  section  describes  the  basic  mathematical  formulations  of  the  thermal-mechanical 
deformation  process  during  welding,  focusing  on  the  formulations  of  the  stress-strain  constitutive 
relation  as  function  of  temperature  and  the  phase  transformation  processes.  These  formulations 
are  essential  to  the  sensitivity  study  of  the  welding  process  effects.  The  mathematical 
formulations  and  numerical  implementation  of  other  aspects  of  thermal-mechanical-metallurgical 
simulation  of  welding  process  can  be  found  in  the  literature  [for  example.  Refs.  8-103.  The 
treatments  of  the  sensitivity  and  fuzzy-set  analysis  can  also  be  found  in  Refs  14  to  16. 

Uncoupled  thermomechanical  analysis  is  performed.  The  temperature  field  is  assumed  to  be 
independent  of  stress  and  strains.  The  governing  equations  for  the  thermal  and  mechanical 
analysis  are  summarized  subsequently. 


2.1  Thermal  Analysis 

2.1.1  Governing  equations 

The  heat  flow  during  arc  welding  was  modeled  as  a  heat  condition  formulation.  The 
convective  heat  transfer  inside  the  weld  pool  due  to  circulation  of  the  molten  metal  was  not 


7 


explicitly  simulated,  but  the  convective  effect  was  accounted  for  in  the  heat  conduction  model  by 
means  of  effective  thermal  properties  such  as  effective  thermal  conductivity.  The  governing 
equation  for  transient  heat  conduction  analysis  is  given  by: 

V  •  (*V7Xr,  /))  +  Q  =  pCp  ( 1 ) 

where  p  is  the  density  of  the  material,  Cp  the  specific  heat,  k  the  thermal  conductivity,  Q  the 
internal  heat  generation  rate,  and  T  the  temperature,  r  is  the  coordinate  of  the  material,  t  is  the 
time. 


2.1.2  Boundary  conditions 

The  heat  from  the  welding  arc  was  treated  as  a  volumetric  heat  source,  taking  the 
ellipsoidal  distribution  proposed  by  Goldak  et  al.  [11]. 


 6>/3 rjEI 


xjxabc 


(2) 


where  q  is  the  energy  density  (W/mm3)  from  the  weld  arc,  77  is  the  arc  efficiency.  For  the 
FCAW  process  used  to  weld  the  gussets,  an  arc  efficiency  value  of  0.8  was  used.  E  and  /  are  the 
welding  voltage  and  current,  respectively.  The  parameters  a,  b,  c  are  related  to  characteristics  of 
the  electrode,  and  represent  the  arc  energy  distribution  profile  of  the  FCAW  process.  The  origin 
of  the  (x,  y,  z)  coordinate  system  is  fixed  at  the  center  of  the  weld  arc. 

For  the  cross-section  model  used  in  this  study,  the  welding  heat  is  effectively  applied 
simultaneously  over  the  weld  length.  This  necessitates  the  following  transformation  of  the 
Goldak  model: 


6^3rjEI  '{ 
^  ~  H'fnabc  ‘ 


f  (jt-jto)1 


(g-gp)2 

b1 


(yy-io))2 


(3) 


where  t  is  the  heating  time,  and  V  is  the  arc  travel  speed.  The  heat  input  is  expressed  in  the  local 
coordinate  system  (x,  £)  associated  with  individual  weld  pass  as  shown  in  Figure  2.  The  center 
of  the  welding  arc  is  located  at  (i0,  £,  0).  tQ  is  selected  such  that,  at  t=0,  the  arc  exerts  negligible 
heating  at  the  modeling  plane.  The  heating  then  increases  as  the  analysis  proceeds,  and  reaches 
the  peak  value  at  t=ta.  The  heat  input  in  Equation  (2)  was  programmed  via  the  DFLUX  user 
subroutine  interface  in  ABAQUS. 

For  all  exposed  surface  of  the  model,  both  convection  heat  loss  and  radiation  heat  loss 
were  considered: 


q,oss=-h(T-TJ-<je{T*-T:)  (4) 

where  h  is  convection  heat  loss  coefficient,  8  the  radiative  emissivity,  and  a  the  Stefan-Boltzman 
constant. 
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2.1.3  Solidification  model 


Under  the  welding  conditions  characteristic  to  this  study,  the  nonequilibrium 
solidification  process  results  in  formation  of  dendrite  microstructure  in  the  weld  pool 
(constitutional  super-cooling)  [12].  The  release  of  the  latent  heat  of  fusion  during  melting  and 
solidification  can  greatly  affect  the  shape  and  size,  and  thermal  stresses  around  the  weld  pool. 
Since  the  present  work  is  concerned  with  the  temperature  around  the  weld  pool,  it  is  important  to 
have  a  more  realistic  model  for  the  release  rate  of  latent  heat  in  order  to  better  simulate  the 
temperature  distribution.  In  this  study,  the  release  of  latent  heat  of  fusion  is  related  to  the 
microscopic  kinetics  of  the  solidification  process  in  the  weld  pool. 

The  release  of  latent  heat  during  solidification  is  treated  as  an  internal  heat  generation 
term  in  the  heat  flow  in  the  governing  equation  for  heat  transfer: 


k  for  Ts<T<Tl 
[0  for  Ts  >T,  or  T  >TL 


(5) 


where  q  is  the  release  rate  of  latent  heat  of  fusion,  and  Ts  and  T*.  are  the  solidus  and  liquidus 
temperature  of  the  alloy. 

The  rate  of  latent  heat  release  is  assumed  to  be  proportional  to  the  rate  of  solid  fraction 
change  in  the  solidification  temperature  range  and  can  be  written  as: 


r  df, 

q>-L«  (6) 

where  L  is  the  volumetric  latent  heat  of  fusion  (J/mm3). 

For  substitutional  alloying  elements,  the  Scheil  equation  can  be  used  to  describe  the 
fraction  of  solid  as  a  function  of  the  solid  composition  at  the  solid-liquid  interface: 

c;=*c0(i  -/,)«*-»  (7) 

DuPont  [13]  developed  a  pseudo  Fe-Mo  binary  phase  diagram  for  AL-6XN,  in  which  the 
solid  interface  composition  is  related  to  the  temperature  at  the  solid-liquid  interface.  Taking 
advantage  of  this  relation,  the  rate  of  latent  heat  release  can  be  written  as: 


=r^L=r^L<dT 

q‘  dt  dC\  dT  dt 

«.(!-*)  L 


dJ_ 

dt 


All  the  material  parameters  in  the  above  equation  are  given  in  Table  1. 
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2.2  Mechanical  Analysis 


The  welding  thermal  cycle  causes  metallurgical  transformations  to  take  place  in  many  high- 
performance  engineering  alloys  of  interest  to  the  Navy.  Those  transformations  occur  on  the 
microstructural  level  (grain  length  scale)  and  manifest  themselves  macroscopically  in  the  form  of 
latent  heat  generation  as  well  as  additional  dilatational  and  deviatoric  deformations.  The  latent 
heat  affects  temperature  history  during  welding  and,  hence,  indirectly  affects  the  residual  stresses. 

Therefore,  in  addition  to  the  elastic,  plastic  and  thermal  strains  in  a  classical  thermomechanical 
deformation  problem,  the  solid-state  phase  transformations  during  welding  give  rise  to  two 
additional  strains.  The  total  strain  rate  can,  therefore,  be  decomposed  into  five  components  as 
follows  (Ref.  17): 


&=  &  +  tSP  +dih  +&  +$p 


(8) 


The  components  on  the  right-hand  side  of  Eq.  (8)  correspond  to  elastic,  plastic,  thermal, 
transformation,  and  transformation  plasticity  strain  rates,  respectively.  In  Eq.  (8),  as  well  as  in 
the  sequel,  boldface  letters  and  symbols  refer  to  second  order  tensors.  It  has  been  shown  [see,  for 
instance.  Refs.  18,  19]  that  the  strain  components  e*  and  s1’’  have  a  major  effect  on  the  welding 
residual  stress  fields  and,  therefore,  should  be  accounted  for  in  any  reliable  analysis  of  those 
fields. 

The  elastic  strain  is  modeled  using  the  isotropic  Hooke's  law  with  temperature-dependent 
Young's  modulus  and  Poisson's  ratio.  For  the  plastic  strain  component  a  rate-independent  plastic 
model  is  employed  with  the  following  features:  kinematic  hardening,  the  von  Mises  yield  surface, 
and  temperature-dependent  material  properties.  Kinematic  hardening  is  considered  a  significant 
feature  since  material  points  typically  undergo  both  loading  and  unloading  during  the  weld 
thermal  cycle  [Ref.  20].  Denoting  the  back  stress  tensor  by  a  and  the  yield  function  by  <}>,  the 
plastic  strain  can  be  expressed  in  the  form: 


afr(g-a) 


(9) 


where 


(10) 


In  Eq.l  1,  is  the  equivalent  plastic  strain  rate: 


(13) 
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In  Eq.  (10),  cto  is  the  initial  yield  stress  surface  which  depends  on  both  temperature  and 
microstructure  constituents,  (  )•(  )  indicates  inner  product  of  two  tensors,  a' is  the  deviatoric 
stress  tensor: 


a'  =  CT-ohI; 


1 

oh  =-tr  a 
3 


(14) 


where  Oh  is  the  hydrostatic  stress,  I  is  the  second-order  identity  tensor,  and  tr  denotes  the  trace  of 
a  tensor.  Similarly,  in  Eq.  (10),  a'  is  the  deviatoric  part  of  a.  At  high  temperatures  the  stress- 
strain  curves  typically  exhibit  pronounced  curvatures  and,  therefore,  a  nonlinear  kinematic 
hardening  rule  is  used  to  describe  the  evolution  of  a.  The  rule  is  based  on  the  work  of  Lemaitre 
and  Chaboche  [Ref.  14]  and  reads: 


+  (15) 

CT0  H  dT  v  ' 

where  H  is  a  temperature  and  microstructure  dependent  hardening  coefficient  and  y  is  a  recovery 
coefficient. 

The  thermal  strain  rate  component  is  given  by  the  relation: 

&  =  (16) 

where  ax  is  the  coefficient  of  thermal  expansion  (CTE)  which  depends  on  temperature  as  well  as 
the  different  phases  present.  The  latter  dependence  is  modeled  using  the  linear  rule  of  mixtures. 
The  CTE  of  the  liquid  phase  is  taken  to  be  zero  to  suppress  the  thermal  expansion  at 
temperatures  above  the  melting  point. 

The  transformation  strain  str  is  purely  dilatational  and  results  from  contraction  and  expansion  of 
the  material  during  allotropic  phase  transformations  (such  as  formation  of  the  austenite  on 
heating  and  decomposition  of  the  austenite  on  cooling  of  ferritic  alloy  steels).  This  strain 
component  can  be  expressed  mathematically  by  the  following  simple  empirical  form  [Ref.  17]: 

(17) 

In  Eq.  (17),  rj  is  a  material  parameter  whose  value  can  be  determined  from  the  total  material 
contraction  and  expansion  induced  by  the  transformation  and  obtained  from  dilatometer 
measurements.  $  is  the  transformation  rate  (l&=  ^Aor  &)• 

The  strain  component  str  is  deviatoric  and  is  linearly  proportional  to  the  stress.  It  occurs  only 
during  the  period  of  martensitic  transformation,  and  is  irrecoverable  upon  removal  of  the  stress. 
This  component  can  be  expressed  using  the  following  empirical  relation  (Ref.). 
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(18) 


«£p=|Mi-4m)^  o' 


where  X  is  a  material  parameter. 

As  shown  above,  the  effects  of  metallurgical  transformations  manifest  in  two  ways  -  the 
additional  two  strain  components  in  Eq.  (8),  and  the  dependency  of  the  yield  stress  surface  on  the 
microstructure  constituents  present  at  the  material  point  in  Eq.  (10)  and  (15). 

3.  Finite  Element  Discretization 

The  heat  transfer  equations  emanating  from  a  finite-element  semi-discretization  are 
integrated  using  an  implicit  backward  difference  scheme  to  generate  the  time  history  of  the 
temperatures.  The  mechanical  response  during  welding  is  them  calculated  by  using  a  rate- 
independent,  small  deformation  thermo-elasto-plastic  material  molded  with  temperature- 
dependent  material  properties. 


4.  Sensitivity  Analysis 

The  sensitivity  coefficients,  which  are  the  derivatives  of  the  various  thermal  and 
mechanical  response  quantities  with  respect  to  the  material  parameters,  are  evaluated  buy  using 
the  direct  differentiation  method  in  conjunction  with  the  automatic  differentiation  software 
facility  ADIFOR.  The  sensitivity  coefficients  obtained  by  ADIFOR  were  validated  by 
comparing  them  with  those  obtained  by  finite  difference  approximations,  using  different  finite 
difference  intervals. 


5.  Numerical  Studies 

A  number  of  studies  have  been  reported  for  the  temperature,  residual  stress-time  histories 
and  their  sensitivity  coefficients  for  a  double  fillet  conventional  welding  of  a  stiffener  and  a  base 
plate  made  of  stainless  AL-6XN  (Ref.  16).  Herein,  the  sensitivity  of  the  temperature  with 
respect  to  microstructural  parameters  is  presented. 

5.1  Material,  Geometry  and  Welding  conditions 

Alloy  AL-6XN,  a  super  austenitic  stainless  steel,  was  used  in  this  study.  It  was  chosen 
because  microsegregation  and  solidification  cracking  have  been  two  major  concerns  in  welding 
AL-6XN  [12].  Table  2  shows  the  nominal  composition  of  AL-6XN  used  in  the  study.  It  is  noted 
that  the  alloy  consists  of  mostly  substitution  alloying  elements,  thereby  the  Scheil  equation, 
would  be  adequate  to  treat  the  micro-segregation  and  release  of  the  latent  heat  of  fusion  during 
solidification. 

The  weld  joint  geometry  used  in  this  study  is  a  double  fillet  Tee-Joint,  representing  a 
stiffened  structure  panel  commonly  used  for  Navy  ship  structure  constructions.  The  dimensions 
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of  the  weld  joint  are  given  in  Figure  2.  The  bottom  panel  plate  was  20-in  long,  20-in  wide, 
whereas  the  vertical  stiffener  was  8-in  high.  Both  the  panel  plate  and  the  stiffener  were  '/2-in 
thick.  The  selection  of  this  particular  panel  configuration  was  consistent  with  previous  studies 
sponsored  by  US  Navy  on  ship  panel  structures  [21,22].  All  welding  conditions  used  in  this 
study  were  also  based  on  the  previous  experimental  welding  tests.  They  are  summarized  below. 

The  stiffened  thin-section  panel  was  welded  using  flux  cored  arc  welding  (FCAW) 
process.  The  double-fillet  welds  were  made  with  twin-torch  configuration:  two  weld  torches 
were  used  to  simultaneously  deposit  the  two  weld  beads  on  two  sides  of  the  stiffener.  The  twin- 
touch  arrangement  is  commonly  used  in  the  shipyard  to  increase  the  productivity  and  to  maintain 
symmetry  about  the  stiffener  for  distortion  control.  The  baseline  welding  parameters  included 
the  followings:  welding  voltage  -  24V,  welding  current  -  245 A,  and  welding  travel  speed  -  14 
in/min.  The  two  welding  guns  were  arranged  such  that  the  second  welding  gun  was  trailing  the 
first  one  by  3.5  in  during  welding.  The  nominal  fillet  size  was  !4  in  wide  and  'A-in  high.  The 
stiffener  was  tack-welded  to  the  panel  at  both  ends  before  the  fillet  welds  were  made. 

5.2  Finite  element  model 

The  two-dimensional  cross-section  model  was  used  for  studying  the  effect  of  alloy 
composition  variations  on  the  temperature  field  during  the  weld  metal  solidification.  In  order  to 
appropriately  simulate  the  release  of  latent  heat  during  solidification  based  on  the  solidification 
kinetics  discussed  in  Section,  the  element  size  in  the  weld  region  needs  to  be  in  the  range  of  0.1 
to  0.2  mm.  A  full  3-dimensional  model  with  such  small  elements  would  be  too  large  to  solve 
practically  with  the  computational  power  available  to  the  program. 

Because  of  its  computational  efficiency,  a  two-dimensional  cross-section  model  was  used. 
This  model  utilized  only  a  cross-sectional  plane,  perpendicular  to  the  weld  length,  of  the  T-joint 
structure  in  the  analysis.  The  appropriate  use  of  the  cross-sectional  model  for  long,  straight 
welds  such  as  those  of  this  study  has  been  well  documented  in  many  studies  in  the  past.  The 
cross-sectional  model  is  generally  accurate  when  the  quasi-steady  state  is  established  in  the 
middle  portion  of  the  welded  structure,  thereby  adequate  for  comparative  evaluations  of  the  alloy 
composition  variations  in  this  study. 

The  2-dimensional  cross-section  model  consists  of  approximately  8600  element.  The  T- 
joint  was  discretized  with  4-noded  quadrilateral  linear  isoparametric  elements.  The  linear 
element  is  more  suitable  for  handling  the  latent  heat  effects  associated  with  welding 
solidification  [23].  Figure  4  shows  the  finite  element  mesh  around  the  double-fillet  welds  of  the 
2-dimensional  model. 

For  evaluation  of  the  sensitivity  of  mechanical  properties  and  welding  process  parameters, 
a  full  3-dimensional  model  was  used.  The  T-joint  was  discretized  with  8-noded  linear 
isoparametric  elements.  The  element  size  in  the  weld  region  was  about  1-mm.  The  entire 
stiffened  panel  consisted  of  approximately  74,000  elements  and  90,000  nodes. 

The  calculations  were  carried  out  with  the  general-purpose  code  ABAQUS  (Version  6.3) 
[24],  enhanced  with  a  set  of  proprietary  user  subroutines  to  deal  with  the  special  issues 
associated  with  welding  process  simulations  that  have  been  described  in  the  previous  sections. 
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5.3  Selection  of  material  parameters  and  their  variation  ranges 

This  study  focused  on  the  sensitivity  study  of  the  material  parameters  mostly  affecting 
the  solidification  kinetics  during  welding.  The  parameters  consisted  of  the  thermophysical 
properties  and  the  solidification  parameters  as  defined  in  Table  .  Their  effects  on  the 
temperature  near  the  solidification  temperature  range  were  analyzed. 

The  thermophysical  properties  (thermal  conductivity,  specific  heat)  were  temperature 
dependent  and  were  taken  from  a  recent  study  by  Nied  et  al  [25].  They  are  provided  in  Figure  6. 
The  thermal  conductivity  above  the  melting  temperature  is  artificially  increased  to  compensate 
for  the  convection  heat  transfer  effect  in  the  molten  weld  pool. 

For  the  sensitivity  analysis,  variations  of  the  above  thermophysical  properties  at  1600K 
(near  the  nominal  solidus  temperature  of  AL-6XN)  was  made  at  ±0.1%,  0.5%  and  1%  to  study 
their  effects  on  the  temperature  field  in  the  solidification  temperature  range.  The  solidification 
parameters  in  Table  were  also  symmetrically  varied  at  ±0.01%,  0.05%,  0.1%,  and  0.5%. 

6.  Results  and  Discussions 


After  the  welding  heat  flow  analyses  were  performed  using  the  finite  element  model,  the 
thermal  histories  at  eight  points  near  the  weld  joint  were  chosen  to  perform  the  sensitivity 
analysis  to  assess  the  effects  of  various  material  parameters  on  the  transient  temperature 
distribution  around  the  molten  weld  pool.  These  eight  points  are  identified  in  Figure  7.  Figure  7 
shows  the  typical  thermal  histories  at  four  points  on  the  right  weld  (Points  5  to  8  in  Figure  6). 
There  are  noticeable  second,  lower  peak  in  the  plots,  associated  with  the  heating  from  the  trailing 
welding  arc  on  the  other  side  of  the  web  plate.  Also  noticeable  in  the  plot  was  the  drastic  change 
in  the  cooling  rate  at  the  center  point  of  the  weld  surface  (Point  6)  due  to  the  release  of  the  latent 
heat  as  the  weld  is  solidified. 

The  finite  difference  sensitivity  analysis  showed  very  good  convergence,  with  decreasing 
the  finite  difference  interval,  for  all  the  material  parameters  studied  here.  As  an  example.  Figure 
9  compares  the  calculated  normalized  sensitivity  of  temperature  with  respect  to  the  nominal  alloy 
composition  for  four  different  finite  difference  intervals  (0.01%,  0.05%,  0.1%,  and  0.5%)  at  the 
center  point  of  weld  surface  (Point  6  in  Figure  7).  Very  good  convergence  was  obtained  for  the 
entire  welding  time  considered.  It  is  also  important  to  notice  that  the  maximum  sensitivity 
occurred  when  the  temperature  was  high,  suggesting  that  the  variations  in  solidification 
parameters  has  the  greatest  influence  in  the  solidification  temperature  range  that  is  critical  to  the 
formation  of  weld  defects. 


7.  Proposed  Future  Work 

Based  on  the  present  study,  future  work  should  focus  on  investigating  the  effects  of 
welding  process  conditions  on  the  microstructural  variations  and  residual  stresses  of  welded 
structures,  by  utilizing  the  computational  welding  simulation  environment  developed  in  the 
present  ONR  study.  The  porposed  research  and  development  activities  include  the  following 
major  efforts: 
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(1)  Refine  the  microstructural  models  to  cover  both  diffusive  and  non-diffusive  phase 
transformation  mechanisms,  and  integrate  these  microstructural  models  into  the 
computational  welding  simulation  environment. 

(2)  Establish  an  effective  computational  methodology  for  sensitivity  analysis  of  the  residual 
stresses  to  variations  of  various  welding  process  conditions.  Well-established  methods 
in  probability  theory  can  be  used  to  characterize  uncertainty  in  process  conditions. 
However,  for  situations  where  the  uncertainty  is  of  a  non-stochastic  nature  (e.g., 
uncertainty  arising  from  imprecision  or  fuzziness),  the  fuzzy  set  approach  has  been 
suggested  to  be  more  appropriate  for  modeling  uncertainty. 

(3)  Study  the  sensitivity  of  key  welding  process  parameters  in  several  promising  welding 
techniques  for  residual  stress  mitigation  such  as  active  quenching,  low  temperature 
phase  transformation  filler  metal,  and  low-stress  and  no-distortion  welding.  Both  arc 
welding  and  laser  welding  will  be  studied. 

(4)  Incorporate  the  local  microstructure  and  property  gradients  in  the  HAZ  and  weld  region 
into  dynamic  failure  and  damage  analysis  of  metallic  structures. 

The  proposed  uncertainty  analysis  will  provide  upper  and  lower  bounds  (i.e.,  envelopes)  for  the 
welding  residual  stresses  and  distortions  as  functions  of  the  variability  in  welding  process 
conditions.  These  bounds  can  help  achieve  the  following: 

•  Better  understanding  of  the  mechanical  response,  fracture  behavior,  and  failure  of  welded 
structural  components; 

•  Better  fitness-for-purpose  assessments  of  welded  structural  components;  and 

•  Optimization  of  the  weld  process  selection  and  control  to  reduce  the  residual  stresses. 
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Table  1  Material  parameters  used  in  this  study  (After  DuPont  |13|) 


Parameters 

Symbol,  unit 

Value 

Nominal  Alloy  Composition 

C0,  Wt-% 

6.3 

Maximum  solid  solubility 

CSE ,  Wt-% 

9.4 

Solute  partition  ratio 

(dendrite  core  /  nominal  alloy 

composition) 

K 

0.66 

16 


Solidus  slope 

ms ,  K/Wt-% 

-17 

Melting  Temperature 

To,0  C 

1490 

Density 

p,  g/mm3 

8.06  x  10'3 

Volumetric  Latent  Heat  of 

Fusion 

L,  J/mm3 

2.1 

Table  2  Nominal  composition  of  AL-6XN,  in  wt-%. 


Fe 

Ni 

Cr 

Mo 

Other 

47 

24.4 

20.9 

6.3 

0.9 
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Figure  2  Local  coordinate  system  used  to  define  the  welding  heat  input. 


Figure  3  Dimensions  of  the  T-joint  used  in  this  study 
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Figure  4  The  twin-welding  torch  arrangement  used  in  making  the  two  fillet  welds 
on  the  T-joint.  The  welding  torch  on  the  left  side  follows  the  right  side 
torch  by  3.5  in. 


Figure  5  Finite  element  mesh  near  the  weld  region  of  the  2-dimensional  cross- 
section  model 
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Figure  6  Temperature-dependent  thermal  conductivity,  specific  heat  and 
convective  heat  loss  coefficient  used  in  this  study  (Ref.  16) 
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Figure  7  Selected  eight  locations  around  the  two  fillet  welds  for  quantifying  the 
sensitivity  of  material  parameters  on  the  temperature  variation  during  welding 
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Figure  8  Variation  of  temperature  as  function  of  welding  time  at  four  different 
locations  (Points  5  to  8)  shown  in  Figure  7. 


time  (secs) 

Figure  9  Time  history  of  the  “normalized”  sensitivity  of  temperature  with  respect 
to  the  nominal  alloy  composition  at  the  center  of  weld  surface  (Point  6  of 
Figure  7) 
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Figure  10  Time  histories  of  the  first-order  sensitivity  coefficient  of  the  temperature 
with  respect  to  the  solidification  parameters  as  the  center  point  of  weld 
surface  (Point  6  of  Figure  7). 
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Appendix  II 

Abstracts  of  Publications 


1 .)  Sensitivity  Analysis  of  the  Thermomechanical  Response  of  Welded  Joints 

A  computational  procedure  is  presented  for  evaluating  the  sensitivity  coeffcients 
of  the  thermomechanical  response  of  welded  structures.  Uncoupled  thermo-mechanical 
analysis,  with  transient  thermal  analysis  and  quasi-static  mechanical  analysis,  is 
performed.  A  rate  independent,  small  deformation  thermo-elasto-plastic  material  model 
with  temperature  dependent  material  properties  is  adopted  in  the  study.  The  temperature 
field  is  assumed  to  be  independent  of  the  stresses  and  strains.  The  heat  transfer  equations 
emanating  from  a  finite-element  semi-discretization  are  integrated  using  an  implicit 
backward-difference  scheme  to  generate  the  time  history  of  the  temperatures.  The 
mechanical  response  during  welding  is  then  calculated  by  solving  a  generalized  plane 
strain  problem.  First-  and  second-order  sensitivity  coeffcients  of  the  thermal  and 
mechanical  response  quantities  (derivatives  with  respect  to  various  thermo-mechanical 
parameters)  are  evaluated  using  a  direct  differentiation  approach  in  conjunction  with  an 
automatic  differentiation  software  facility.  Numerical  results  are  presented  for  a  double 
fillet  conventional  welding  of  a  stiffener  and  a  base  plate  made  of  stainless  steel  AL-6XN 
material.  Time  histories  of  the  response  and  sensitivity  coeffcients,  and  their  spatial 
distributions  at  selected  times  are  presented. 


2.)  Uncertainty  Analysis  of  Welding  Residual  Stress  Fields 

A  fuzzy-set  approach  is  applied  in  conjunction  with  a  finite  element  procedure  to 
study  the  effect  of  uncertainty  in  material  properties  on  welding  residual  stress  fields.  The 
study  focuses  on  steels  that  undergo  Martensitic  transformations  during  the  cool-down 
part  of  the  weld  thermal  cycle.  Uncoupled  thermo-mechanical  analysis  is  performed. 

Filler  metal  deposition,  various  convective-radiative  energy  losses,  and  latent  heat 
generation  are  included  in  modeling  of  the  welding  process.  A  rate-independent  thermo- 
elasto-plastic  material  model  with  nonlinear  kinematic  hardening  is  adopted  in  the  study. 
Account  is  taken  of  the  metallurgical  transformations  that  occur  during  the  entire  weld 
thermal  cycle.  The  material  model  is  validated  by  comparing  its  predictions  with  existing 
experimental  measurements  of  the  welding  residual  mechanical  fields.  Fuzzy  logic  is  then 
used  to  study  the  variability  in  the  residual  stress  field  due  to  pre-selected  ranges  of 
variation  in  material  parameters  associated  with  the  Martensitic  transformations.  It  is 
found  that  the  variations  in  material  parameters  considered  have  a  significant  effect  on 
the  welding  residual  stress  field. 
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3.)  A  Fuzzy-Set  Analysis  for  a  Dynamic  Thermo-elasto-viscoplastic  Damage 
Response 

A  computational  fuzzy-set  based  approach  is  used  to  study  variability  (uncertainty) 
in  the  mechanical  response  of  dynamically-loaded  simple  structural  components.  The 
mechanical  response  is  described  by  a  thermo-elasto-viscoplastic  material  model  that 
accounts  for  damage.  The  uncertainty  in  the  response  is  associated  with  pre-selected 
variations  in  material  parameters.  Sensitivity  studies  are  performed  to  identify  the 
material  parameters  which  have  the  most  significant  effect  on  the  response.  Those 
parameters  are  then  assumed  to  be  fuzzy,  and  time  histories  of  the  possibility 
distributions  of  the  response  quantities  are  constructed.  Numerical  results  are  presented 
for  a  titanium  alloy  Ti-6A1-4V. 


4.)  Sensitivity  Analysis  for  Failure  and  Damage  in  Dynamically  Loaded  Tensile 

Bars. 


A  computational  procedure  is  presented  for  evaluating  the  sensitivity  coefficients 
of  porous  viscoplastic  solids  under  dynamic  loading  conditions.  The  effects  of  finite 
strains,  material  strain  and  strain-rate  hardening,  and  the  thermal  softening  due  to 
adiabatic  heating  are  incorporated  into  the  formulation.  A  critical  void  volume  fraction 
criterion  is  used  to  identify  the  initiation  of  failure.  The  equations  of  motion  emanating 
from  a  finite  element  semi-discretization  are  integrated  using  an  explicit  central 
difference  scheme.  First-  and  second-order  sensitivity  coefficients  of  the  response 
quantities  (derivatives  with  respect  to  various  material  parameters)  are  evaluated  using  a 
direct  differentiation  approach  in  conjunction  with  an  automatic  differentiation  software 
facility.  Numerical  results  are  presented  for  tensile  specimens  subjected  to  impact  loading. 
Both  axisymmetric  and  plane  strain  states  are  considered.  The  first-  and  second-order 
sensitivity  coefficients  are  generated  by  evaluating  the  derivatives  of  the  response 
quantities  with  respect  to  various  macroscopic  and  microscopic  material  parameters. 

Time  histories  of  the  response  and  sensitivity  coefficients,  and  their  spatial  distributions 
at  selected  times  are  presented. 
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Appendix  III 


Screenshots  of  Presentation  on  Residual  Stresses  and  Distortions 
in  Large  Welded  Structures 
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Future  Work 

Assess  the  reliability  of  the  materials  models 
incorporating  microscopic  solidification  and 
metallurgical  transformation  (with  Oak  Ridge, 
Lehigh  University  and  Penn  State  teams) 

-  Systematic  verification,  validation  and 
calibration  of  models 

„ j  Explore  processing  methods  for  reducing 
residual  stresses  and  distortions  in  stainless 
steel  ADH  fabrication  (with  Oak  Ridge,  Lehigh 
University,  Penn  State,  and  Navy  Lab  teams) 

-  Effective  use  of  sensitivity  information 


Mam  Menu  I 

Hi—  -  1 

Future  Work 

Apply  the  methodologies  developed  to  large 

welded  structures  (with  Oak  Ridge,  Navy  Lab 

and  other  teams)  _ 

‘  A?  k 

. .  •  i! 
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Outline 

|  Oojecriyes  and  Scope 

Sensitivity  Information 
Mathematical  Modeling 
Computational  Procedure 
Mumerical  Results 
Summary 


Overall  goals 


Develop  models  and  establish  effective  method¬ 
ologies  for  the  reliable  prediction  of  residual 
stresses,  distortions  and  microstructure  changes 
in  welded  components 


|  Study  effect  of  residual  stresses  and  micro- 
structural  changes  on  both  brittle  and  ductile 
fracture  behavior  * 
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Objectives 

j  Assess  sensitivity  of  thermomechanical  response 
and  damage  of  welded  components  to  various 
microstructural,  and  other  material,  geometric 
and  welding  parameters 

j  Study  effect  of  variability  in  material  and  other 
parameters  on  welding  residual  stresses  and 
distortions 


Importance 

Better  understanding  of  mechanical  response, 
fracture  behavior  and  failure  of  welded 
structural  components 

.j  Better  fitness-  for-  purpose  assessments  of  welded 
components 


yj  Optimization  of  the  weld-  process  selection  and 
control  to  reduce  residual  stresses  and  distortions 
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Scope 

jl  Stainless  steel :  AL-6XN 
^  Heat  source  model 


^  Heat  transfer  model  incorporating  microstructure 
solidification  (Z.  Feng) 


Rate  independent  thcrmo-elasto-pl3Stic 
constitutive  model 


Uses 

Assess  effect  of  uncertainties  | 

^Predict  changes  in  respoi 


nse 


y^Help  in  identifying  regions  of 
the  structure  requiring 
models  with  different 
length  scales 

^ identify  changes  in  para¬ 
meters  required  to  achieve 
certain  performance 
characteristics 


*  Material  and  geometric 
parameters 


•Weld  characteristics 


Microstructural 
parameters 


\jDctcn 

optimi 


PW  complete  j 
M»*n  Mrrmj 


mine  sear 
lization  te< 


search  direction 
techniques 


1 
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Approaches  for  Sensitivity  Analysis 

Direct  Differentiation 

-  Automatic  Differentiation  ( eg ADIFORJ 
-Jump  discontinvities  can  occur. 

Regularization  poses  a  challenge 

yj  Finite  Differences 

-  F.D.  increment  must  be  chosen  to  avoid 
both  truncation  and  roundoff  errors 


^  Mixed  Hybrid  Approaches 


Monlinear  Response: 


(  l  (A\  A,  t)  |  0 


)  X  !  -  fundamental  unknowns  (response  vector) 
/.  *  typical  system  parameter 


time 

Sensitivity  Coefficients: 
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T  riHomition 
caused  try  stress 


Thermal 


Temperature 

field 


Metallurgical 


Microstructura! 
state  field 


Thermal  Mechanical  and  Metallurgical 


MunMmi 
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Pseudo-binary  phase  diagram  for  AL-6XN 


Thermal  6  0  s 
Heal  Source  Data 
Thermomectufiica!  Prop 


Alloy  Composition 
Processing  Parameter* 


Component 
geometry  / 
Discrete  Model 


rinalP 

'£221 


Transformation 


Mechanical  B.C.'tf 
Loading  Data 
Mechanical  Properties 


e  g..  Liquidus  Slopes 
and  phase  boundaries 


Integrated  Approach  for  Relating  Weld  Performance  to 


*  — ~  PHI - 

il'aal 

Performance  Model 


Two-stepsimulation 
-  During  welding  -  Process  model 
-After  welding  -  Performance  model 

Process  model  should  adequately  cover  the 
thermal  -  mechanical  -  metallurgical  aspects 
in  welding 

Spatial  distributions  of  thermomechanical 
properties  are  required  as  input  to  performance 

niraali 
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Welding  Heat  Input 
Double  Ellipsoid  Model  (Goldak  et  a I.) 


q  -  power  density  in  front  (and  rear)  ellipsoids 


a  *  b,  f  «  0.6  before  the  torch  passes  the  analysis  plane 
a  -  4b.  f  -  1 .4  after  the  torch  passes  the  analysis  plane 
b,  c  weld  width  penetration 
f  =»  fractions  of  heat  deposited 
q.  U,  I  ~  welding  efficiency,  voltage  and  amperage 
x.  y,  ^  *  local  coordinates  of  double  ellipsoid  model 
v  *  torch  travel  speed 


1  Heat  Transfer  Model 

(Incorporating  Microstructure  Solidification) 

O  AL6XN  (and  other  stainless  steels)  solidifies 
under  unequilibrium  conditions  due  to  the 
rapid  cooling  in  the  weld  pool 

^  Under  the  welding  conditions  characteristic  to 
this  study ,  the  unequilibrium  solidification 
process  results  in  formation  of  dendrite 
microstructure  in  the  weld  pool  (constitutional 
super-cooling)  (DuPont,  private  communications! 


„ j  The  grain  morphology  is  simplified  as  cellular 
cells 
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Heat  Transfer  Model 

(Incorporating  Microstructure  Solidification  -  Feng) 


[  V.(kVT)+Q=pCpT 

k  -  temperature -dependent  conductivity 
The  release  of  latent  heat  from  formation 
of  dendrite  microstructure  is  treated  as  an 
internal  heat  generation  term  in  the  heat 
flow,  Q 

p  -*  mass  density 

Cp  *  temperature -dependent  specific  heat 


L 


Heat  Transfer  Model 

(Incorporating  Microstructure  Solidification  -  Feng) 


0=  </,  for  t<t, 

“  0  fnr  7  J\.  nr  T  >  T {  ® 

The  release  of  latent  heat  during  weld 
solidification  is  determined  from  the  rate 
of  dendrite  formation 

Jclt  =  Rate  of  latent  heat  release 
Ts  =  Solidus  temperature 
Tf  =  Liquidus  temperature 
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Heat  Transfer  Model 

(Incorporating  Microstructure  Solidification  -  Feng) 

:,'r ^ 


L  ”  total  energy  released  during  solidification  per 
unit  volume  (  volumetric  latent  heat) 
f#  =  fraction  of  solid  during  welding  related  to  the 
composition  at  the  solid-liquid  interface  C$  by 
the  local  solute  redistnbution  equation 

The  solidification  rate  is  linked  to  the  cooling 
rate  through  Fleming's  local  solute  distribution 
equation,  and  phase  diagrams 


Heat  Transfer  Model 


(Incorporating  Microstructure  Solidific 

"  (X  +  ?TCl 

T-Ts«- 

'\-T 

mx(l-A ) 

si: 

:t 

n?s  =  Solidus  slope  C  .  =  Maximum  solid  solubility 

Sh 


C0 *  *  nominal  alloy  composition 


_  hqiudus  slope  mL 
solidiLs  slope  mf 

*  equilibnum  distnbution  coefficient  (partition 
ratio) 

71  =  Solidus  temperature 
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Heat  Transfer  Model 

(Incorporating  Microstructure  Solidification  -  Feng) 

■  Surface  heat  loss  from  convection  and  radiation  are 
included  as  boundary  conditions 

q^hlT-T^-ot^-T:) 

h  *  convection  coefficient 
*  radiation  coefficient 


Cross  Section  of  Fillet  Weld 
Details  of  HAZ  (AL-6XN) 
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Modeling  Metal  Deposition 

First  assign  artificially  soft  properties  to  ail  the 
elements  to  be  welded  later 

^Switch  back  to  the  actual  properties  on 

integration  point  basis  once  it  begins  to  solidify 
from  weld  pool 

Switching  is  done  based  on  temperature  to 
avoid  numerical  convergence  problems 

^Unrealistic  plastic  strains  in  the  weld  pool  due 
to  solid  mechanics  formulation  are  also 
eliminated  upon  switching 


^  Element  rebirth  is  done  on  pass,  not  elements 


Computational  Procedure 
Finite  Element  Implementation 


r  '  =  thermo  metallurgical  strain  tensor 
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Double  Fillet  Conventional  Welding 
of  a  stiffener  and  a  base  plate 


Microstructurat  Parameters 
Nominal  Alloy  Composition,  Co  6.3 

Maximum  solid  solubility,  Cm  9.4 

Partition  ratio  (dendrite  core  /  nominal  alloy 


K  0.66 

m.  -17 

To  1490  °C 

8.06  x  10  1  g/mmJ 
L  2.1  J/mm3 


composition), 

Solidus  slope, 

Melting  Temperature, 
Destiny,  p 

Volumetric  Latent  Heat, 


mm 
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Sm 


;  ; ;  ; i * 


Double  Fillet  Conventional  Welding 
of  a  stiffener  and  a  base  plate 


3D 


Finite 


Element 


Model 


Problem 


/  Number  of  nodes:  90267 
Number  of  elements:  74200 

A'X  I, 

<  t  :  -<f9  ::  ■ 

!U  iJrr  ■  ••• 


' 
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Summary 

,jt  Sensitivity  analysis  used  to  study  effect  of 
variability  of  microstructural  and  other  material 
parameters  on  thermomechanical  response  of 
welded  structures 


Effect  of  microstructure  solidification  is 
incorporated  into  heat  transfer  model 

Rate  independent  thermo-elasto-plastic 
model  used 


Summary 


Interrogative  visualization  capability  d 
to  provide  insight  into  the  response  and 
sensitivity  characteristics 


ity  information  is  used  with  a  rapid 
reanalysis  technique 


Dynamic  meshing  /  Reduction  technique 
developed  for  3-D  simulation  of  large  welded 
structures 


III  a  aM 


r 


MjinMenu 
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Future  Work 


vj  Assess  the  reliability  of  the  materials  models 
incorporating  microscopic  solidification  and 
metallurgical  transformation  (with  Oak  Ridge, 
Lehigh  University  and  Penn  State  teams) 

-  Systematic  verification,  validation  and 
calibration  of  models 

Explore  processing  methods  for  reducing 
residual  stresses  and  distortions  in  stainless 
steel  ADH  fabrication  (with  Oak  Ridge,  Lehigh 
University,  Penn  State ,  and  Navy  Lab  teams) 

-  Effective  use  of  sensitivity  information 


MjmMmu  1 

!■!  _  -A 
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Future  Work 


Apply  the  methodologies  developed  to  large 
welded  structures  (with  Oak  Ridge,  Navy  Lab 
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